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ABSTRACT: The seventh Jerry L. Swedlow Memorial Lecture presents a review of 
some of the technical developments, that have occurred during the past 40 years, which 
have led to the merger of fatigue and fracture mechanics concepts. This review is made 
from the viewpoint of "crack propagation." As methods to observe the "fatigue" 
process have improved, the formation of fatigue micro-cracks have been observed earlier 
in life and the measured crack sizes have become smaller. These observations suggest that 
fatigue damage can now be characterized by "crack size." In parallel, the crack-growth 
analysis methods, using stress-intensity factors, have also improved. But the effects of 
material inhomogeneities, crack-fracture mechanisms, and nonlinear behavior must now be 
included in these analyses. The discovery of crack-closure mechanisms, such as plasticity, 
roughness, and oxide/corrosion/fretting product debris, and the use of the effective stress- 
intensity factor range, has provided an engineering tool to predict small- and large-crack- 
growth rate behavior under service loading conditions. These mechanisms have also 
provided a rationale for developing new, damage-tolerant materials. This review suggests 
that small-crack growth behavior should be viewed as typical behavior, whereas large- 
crack threshold behavior should be viewed as the anomaly. Small-crack theory has unified 
"fatigue" and "fracture mechanics" concepts; and has bridged the gap between safe-life 
and durability/damage-tolerance design concepts. 
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In 1965, ASTM Committee E24 on Fracture Testing of Metallic Materials was formed 
to promote the rapid growth of the field of Fracture Mechanics. Committee E09 on Fatigue 
had already been in existence for nearly two decades. With the recent merger of Committees 
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E09 and E24 to form Committee E08 in 1993, it seems appropriate to explore the merging of 
fatigue and fracture mechanics concepts in a forum that honors and remembers Professor Jerry 
L. Swedlow, who incidentally completed his Ph.D. Thesis also in 1965 [7]. 

Since the 1950’s, events in the naval, nuclear, and aircraft industries have fostered the 
development of the field of fracture mechanics. The failure of the Comet transport jet aircraft 
[2] from fatigue cracks gave rise to treatments of crack propagation using notch-root 
parameters and the stress-intensity factor concept of Irwin [3] and Paris et al. [4,5]. Crack 
propagation theories would eventually form the bridge that would link fatigue and fracture 
mechanics concepts. The notch-root local-stress approach hinged upon the Neuber [6] or 
Hardrath-Ohman [7] equations which related local plastic stresses and strains to the elastic 
stress concentration. Later, Hutchinson [3] and Rice [9] noted some similarities between 
Neuber’s elastic-plastic relation for notches and their solutions for elastic-plastic behavior of 
cracks. Using a notch-root parameter, K(nS ne t, for a sharp notch or crack, McEvily and Dig 
[70] correlated fatigue-crack-growth rates in a very similar manner to the current AK-rate 
concept. Years later, this notch-root parameter was shown to be directly related to the stress- 
intensity factor [5]. But the elegance and simplicity of the stress-intensity factor concept 
rapidly developed into the durability and damage tolerance concepts currently used today to 
design fatigue- and fracture-critical components. The next major link between fatigue and 
fracture mechanics was the discovery of fatigue-crack closure by Elber [77]. The crack- 
closure concept put crack-propagation theories on a firm foundation and allowed the 
development of practical life-prediction methods for variable-amplitude and spectrum loading, 
such as experienced by modern-day commercial aircraft. 

In the mid- 1 970’ s, Pearson [72] and Kitagawa [73] showed that short cracks (less than 
about 0.5 mm in length) grew much faster than long cracks when correlated against the stress- 
intensity factor range. During the next two decades, short- or small-crack research formed the 
final link between fatigue and fracture mechanics. These studies, conducted by many world- 
wide organizations [74,75], the AGARD Structures and Materials Panel [76-73], ASTM 
Committees E9 and E24 [79], NASA and the CAE [20] provided experimental databases and 
analysis methods to perform fatigue analyses on notch components using "crack propagation" 
theories. The small-crack theory (treatment of fatigue as the growth of micro-cracks, 1 to 20 
pm in length) has been applied to many engineering materials with reasonable success. 
Although this review will concentrate mainly on a fracture-mechanics viewpoint, the local 
notch-root stresses and strains from classical fatigue analyses are the driving forces behind the 
initiation and growth of small cracks at material discontinuities or manufacturing defects. The 
merging of fatigue and fracture mechanics concepts will provide industries with a unified 
approach to life prediction. Small-crack theory can now be used to assess the influence of 
material defects and manufacturing or service-induced damage on fatigue life behavior. This 
approach will ultimately improve the reliability and economic usefulness of many structures. 

The review will begin with some fatigue and fatigue-crack growth observations that have 
set the stage for the treatment of fatigue from a fracture-mechanics viewpoint. This treatment 
hinges strongly upon whether fracture-mechanics parameters can be used to model micro- or 
small-crack growth rate behavior. The development of the stress-intensity factor and some 
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nonlinear fracture-mechanics parameters, such as the J and T* integrals, and their application to 
small-crack behavior will be discussed. The application of linear elastic fracture mechanics, i.e. 
the stress-intensity-factor range, AK, to the "small or short" crack-growth regime has been 
actively studied and questioned for more than two decades. Various nonlinear crack-tip 
parameters and crack-closure effects were introduced to help explain the differences between 
small- and large-crack growth rate behavior. A key element in these nonlinear crack-tip 
parameters is crack closure. A brief survey of the finite-element and finite-difference analyses 
that have been conducted to study the fatigue-crack growth and closure processes; and some 
typical results will be given. A review of some of the more popular yield-zone models, 
empirical crack-closure models, and the analytical crack-closure (modified Dugdale or strip- 
yield) models will be discussed. The application of some of these models to predict crack 
growth under aircraft spectrum load histories will be presented. Constraint or three- 
dimensional stress state effects play a strong role in the fatigue initiation and crack-growth 
process. For example, plasticity-induced crack closure (yielded material at the crack tip and in 
the wake of the advancing crack) is greatly affected by plane-stress or plane-strain behavior. 

The most common constraint parameters, and their use in fatigue-crack growth relations, will 
be discussed. The evolution of some of the proposed fatigue-crack-growth rate relations will 
then be reviewed. Some observations on the effects of microstructure, environment, and 
loading on fatigue-crack-growth rate behavior will be discussed. These observations are 
important in developing the intrinsic crack-growth-rate relations to calculate crack growth 
under general cyclic loading. The small-crack growth rate data presented by Pearson [72], and 
enlarged upon by Lankford [27], will be presented and discussed. An analysis of the Pearson 
and Lankford small-crack data reveals an important conclusion about the relevance of the 
large-crack thresholds. The prediction of fatigue life, on the basis of crack propagation from 
microstructural features, such as inclusions or voids, will be presented for several materials and 
loading conditions. A design concept using "small-crack theory" will be discussed. 

This review is necessarily limited in scope and will not be able to fully cover the vast 
amount of research that has been conducted over the past 40 years in the fields of fatigue and 
fracture mechanics. Several excellent books and articles on the "merging of fatigue and 
fracture mechanics concepts" have helped set the stage for this paper. The book by Fuchs and 
Stephens [22] presents a brief history on the subject, the book entitled "Fatigue Crack Growth 
— 30 Years of Progress" edited by Smith [ 23 ] gives excellent reviews on a variety of technical 
subjects, and the paper on "History of Fatigue" by Schutz [ 24 ] gives a historical perspective. 
The author request the readers indulgence and forgiveness if some major events have been 
omitted unintentionally, or if reference is not made to all of those who have made significant 
contributions to the subject. As pointed out by Paris in the Third Swedlow Lecture [25], 
"History has a strong tendency to be one man’s personal recollection of important events ...". 
This review is no different. 
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FATIGUE AND FATIGUE CRACK GROWTH OBSERVATIONS 


The fatigue life, as presented by Schijve [26], is divided into several phases: crack 
nucleation, micro-crack growth, macro-crack growth, and failure, as shown in Fig. 1. Crack 
nucleation is associated with cyclic slip and is controlled by the local stress and strain 
concentrations, and notch constraint. Although the slip-band mechanism of crack formation 
may be necessary in pure metals, the presence of inclusions or voids in engineering metals 
will greatly affect the crack-nucleation process. Micro-crack growth, a term now referred to 
as the "small-crack growth" regime, is the growth of cracks from inclusions, voids, or slip 
bands, in the range of 1 to 10 pm in length. Schijve [27] has shown that for polished 
surfaces of pure metals and for commercial alloys, the formation of a small crack to about 
100 pm in size can consume 60 to 80 % of the fatigue life. This is the reason that there is 
so much interest in the growth behavior of small cracks. Macro-crack growth and failure 
are regions where fracture-mechanics parameters have been successful in correlating and in 
predicting fatigue-crack growth and fracture. This review will highlight the advances that have 
been made in the use of the same fracture-mechanics parameters in the treatment of micro- or 
small-crack growth using continuum-mechanics approaches. 
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FIG. 1— Different phases of fatigue life and relevant factors (modified after Schijve, 1979 [27]). 


One of the earliest observations on the mechanism of small-crack growth was made 
by Forsyth [28]. He showed that the initiation and early growth of small cracks can occur 
at a single slip system (Stage I crack growth) in a favorably oriented surface grain, as 
shown in Fig. 2(a). Slip-band cracking is promoted by high stresses and higher alloy 
purity, such as observed in cladding on aluminum alloys. The transition from Stage I to a 
crack-growth mechanism involving multiple slip systems at the crack tip (Stage II) can 
occur at or near the first grain boundary encountered by the crack. As might be expected, 
grain boundaries can have a significant effect on the growth of small cracks. The grain 
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boundaries contribute greatly to the scatter that is observed in small-crack growth rate 
behavior. 

All materials are anisotropic and inhomogeneous when viewed at a sufficiently small 
scale. For example, engineering metals are composed of an aggregate of small grains. 
Inhomogeneities, see Fig. 2(b), exist not only due to the grain structure, but also due to 
the presence of inclusion particles or voids. These inclusion particles are of a different 
chemical composition than the bulk material, such as silicate or alumina inclusions in 
steels. Because of the nonuniform microstructure, local stresses may be concentrated at 
these locations and may cause fatigue cracks to initiate. Crack initiation is primarily a 
surface phenomenon because: (1) local stresses are usually highest at the surface, (2) an 
inclusion particle of the same size has a higher stress concentration at the surface than in 
the interior, (3) the surface is subjected to adverse environmental conditions, and (4) the 
surfaces are susceptible to inadvertent damage. The growth of "natural" surface initiated 
cracks in commercial aluminum alloys has been investigated by Bowles and Schijve [29], 
Morris et al. [30] and Kung and Fine [31]. In some cases, small cracks initiated at 
inclusions and the Stage I period of crack growth was eliminated, as shown in Fig. 2(b). 
This tendency toward inclusion initiation rather than slip-band (Stage I) cracking was 
found to depend on stress level and inclusion content [31]. Similarly, defects (such as tool 
marks, scratches and burrs) from manufacturing and service-induced events will also 
promote initiation and Stage II crack growth, as shown in Fig. 2(c). 

In 1956, Hunter and Fricke [32] conducted rotating beam tests on chemically- 
polished un-notched specimens made of 6061-T6 aluminum alloy. Testing was interrupted 
periodically in order to obtain plastic replicas of the specimen surface. These replicas 
reproduced the surface details and provided direct measurement of cracks. The stress- 





Stages of fatigue process influence of inclusions Influence of service-induced 

or voids or manufacturing defects 

FIG. 2— Stages of the fatigue-crack-growth process (after Forsyth, 1962 [28] ; after Morris, 
Buck and Marcus, 1976 [30]). 
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cycles (S-N) relation between the "first" crack observed and failure is shown in Fig. 3. 
(No crack length was defined for the first crack, but crack length data was presented for 
lengths greater than about 0. 1 mm.) These tests revealed that at high stresses, crack 
propagation was a dominate part of life, whereas at low stresses, near the endurance limit, 
crack nucleation was dominant. Their results on crack-length-against-cycles from 0. 1 to 1 
mm did not show any abnormal behavior (i.e., rates were a monotonically increasing 
function of crack length). 



N, cycles 


FIG. 3— Stress-life curves for rotating beams under constant-amplitude loading (after Hunter 
andFricke, 1956 [ 32 ]). 


Newman and Edwards [77], in an AGARD cooperative test program involving 
several laboratories, found similar results under an aircraft load spectrum, FALSTAFF, on 
chemically-polished notched specimens made of 2024-T3 aluminum alloy. Replicas were 
also used to monitor crack initiation and crack growth from 10 pm to 2.3 mm. The stress- 
life curves for lives to a given crack length are shown in Fig. 4. These results show that 
crack growth is the dominant part of life (about 90 %) for the limited range of stress levels 
tested. Similar results were also found for the Gaussian load spectrum and various 
constant-amplitude loading conditions (R = -2, -1, 0 and 0.5) with stress levels above and 
near the endurance limit. For all loading conditions, 80 to 90 percent of the fatigue life 
was spent as crack growth from a crack length of 20 pm to failure (inclusion-particle 
cluster sizes at the initiation sites ranged from 2 to 7 pm in length). How much of the 
fatigue life is consumed by crack growth from a crack of the inclusion-particle size to 20 
pm is left to conjecture, but it could easily account for 5 to 10 percent of the total life. 
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FIG. 4— Stress-life curves for single-edge-notch tensile (SENT) specimens under aircraft 
spectrum loading (after Newman and Edwards, 1988 [77]). 


Elber [33], in 1968, observed that fatigue-crack surfaces contact with each other 
even during tension-tension cyclic loading. This contact is due to residual plastic 
deformation that is left in the wake of an advancing crack, as illustrated in Fig. 5(a). This 
deformed material contacts during unloading. It is surprising that this observation 
appeared so many years after crack growth was first studied. But this simple observation 
and the explanation of the crack-closure mechanism (or more properly crack-opening) 
began to explain many crack-growth characteristics almost immediately. Since the 
discovery of plasticity-induced closure, several other closure mechanisms, such as 
roughness- and oxide/corrosion/fretting product-induced closure, have been identified. 

The roughness mechanism, discovered by Beevers and his coworkers [34, 55], appears to 
be most prevalent in the near-threshold regime of large-crack growth where the maximum 
plastic-zone sizes are typically less than the grain size [36]. At these low stress levels, 
crack extension is primarily along a single slip system resulting in a Stage I-like mechanism 
and a serrated or zig-zag (± 0 deg.) crack-growth path, as shown in Fig. 5(b). These 
cracks will have mixed-mode (Mode I and II) crack-surface deformations, which provide 
the mechanism for contact between the surfaces during cyclic loading. Cracks growing 
along a non-planar path, such as during overloads in aluminum alloys, will develop surface 
contact and create debris due to fretting and the growth of oxides from the newly created 
crack surfaces, see Fig. 5(c). This debris will cause premature contact, as discussed by 
Paris et al. [57] and Suresh et al. [55]. These new closure mechanisms, and the influence 
of the plastic wake on local crack-tip strain field, have greatly advanced the understanding 
of the fatigue-crack growth process. A brief review of some of the numerical analyses and 
models of the crack-closure phenomenon will be presented later. 
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(a) Plasticity-induced (b) Roughness-induced 
closure closure 


(c) Oxide/corrosion product- 
induced closure 


FIG. 5--Dominant fatigue-crack-closure mechanisms (after Suresh and Ritchie, 1982 [36]). 


STRESS-INTENSITY FACTORS 

The essential feature of fracture mechanics is to characterize the local stress and 
deformation fields in the vicinity of a crack tip. In 1957, Irwin [3, 39] and Williams [40] 
recognized the general applicability of the field equations for cracks in isotropic elastic bodies. 
Under linear elastic conditions, the crack-tip stresses have the form: 

Oij = K (2tu r)' 1/2 fjj(9) + A 2 gjj(e)+ A 3 hij(0)r 1/2 + ... (1) 

where K is the Mode I stress-intensity factor; r and 0 are the radius and polar angle measured 
from the crack tip and crack plane, respectively; A; are constants; fjj(0), gij(0) and hy(0) are 
dimensionless functions of 0. The stress fields for both two- and three-dimensional cracked 
bodies are given by equation (1). After some 30 years, the stress-intensity factors for a large 
number of crack configurations have been generated; and these have been collated into several 
handbooks (see for example Refs. 41 and 42). The use of K is meaningful only when small- 
scale yielding conditions exist. Plasticity and nonlinear effects will be covered in the next 
section. 

Because fatigue-crack initiation is, in general, a surface phenomenon, the stress-intensity 
factors for a surface- or comer-crack in a plate or at a hole, such as those developed by Raju 
and Newman [43, 44], are solutions that are needed to analyze small-crack growth. Some of 
these solutions are used later to predict fatigue-crack growth and fatigue lives for notched 
specimens made of a variety of materials. 
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ELASTIC-PLASTIC OR NONLINEAR CRACK-TIP PARAMETERS 


Analogous to the stress field for a crack in an elastic body, Hutchinson, Rice and 
Rosengren (HRR) [ 8 , 45 ] derived the asymptotic stress and strain field for a stationary crack in 
a nonlinear elastic body. The first term for a power-hardening solid was given by: 


= [E'J/(a o 2 r)] n/(n+1) a o 5j(0,n) 

(2) 

= [E' J/(a 0 2 r)] 1/(n+1) gij(0,n) 

( 3 ) 


where J is the path-independent integral of Rice [ 9 ], E' is the elastic modulus (E 1 = E for plane 
stress or E' = E/(l-v ) for plane strain), g 0 is the flow stress, n is the strain-hardening 
coefficient (n = 0 is perfectly plastic and n = 1 is linear elastic), r and 0 are the radius and polar 
angle measured from the crack tip and crack plane, respectively, and fjj(0,n) and gij(0,n) are 
dimensionless functions depending upon whether plane-stress or plane-strain conditions are 
assumed. 

J and T* Path Integrals 

The J-integral appeared in the works of Eshelby [ 46 ], Sanders [ 47 ], and Cherepanov 
[ 48 ], but Rice [ 9 ] provided the primary contribution toward the application of the path- 
independent integrals to stationary crack problems in nonlinear elastic solids. The J- 
integral, defined in Fig. 6, has come to receive widespread acceptance as an elastic-plastic 
fracture parameter. Landes and Begley [ 49 ] and others have used the J-integral as a 
nonlinear crack-tip parameter to develop crack initiation Ji c values and J-R curves for a 
variety of materials. However, because deformation theory of plasticity, instead of 
incremental theory, was used in its derivation, the J-integral is restricted to limited 
amounts of crack extension in metals. 

In 1967, Cherepanov [ 48 , 50] derived an invariant integral (denoted F) that is valid for 
the case of a moving crack with arbitrary inelastic properties, such as an elastic-plastic material. 
(Both Rice and Atluri have used the symbol T to denote a contour around the crack tip. This 
should not be confused with Cherepanov’s definition of his T-integral.) Atluri and his co- 
workers [57, 52] overcame experimental and numerical difficulties in evaluating the T- 
parameter (denoted as T*, see Fig. 6) for a moving crack. The T*-integral is beginning to -~ 
receive more attention in the literature. 

Cyclic Crack-Tip Parameters 

Both the J- and T*-integrals have been extended to apply to applications involving cyclic 
loading. Dowling and Begley [55] developed an experimental method to measure the cyclic J 
values from the area under the load-against-deflection hysteresis loop that accounted for the 
effects of crack closure. The AJ e gf parameter has been successfully used to correlate fatigue- 
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crack-growth rate data from small- to large-scale yielding conditions for tension and bending 
loads [ 54 \. Similarly, Atluri et al. [ 52 ] have derived the AT* integral (see Fig. 6) for cyclic 
loading and others [55] are beginning to evaluate the parameter under cyclic loading. 


Plastic zone 


Plastic wake 




Rice (1968): Atluri (1982): 

J = ^ (Wn 1 - tj Uj .| )dr T*= J - ^ [W,., - (ay £jj -|)]dV 
r ’ V-V £ 

Dowling and Begley (1976): Atluri, Nishioka and Nakagaki (1984): 

AJ eff = 2 A c / B b AT* = \(AWn 1 - (tj + Atj ) Auj - Atj u it1 ) dr 

r £ 

FIG. 6— Elastic-plastic crack-tip parameters. 


Plastic Stress-Intensity Factors and the Duedale Model 

In 1960, Irwin [ 56 ] developed a simple approach to modify the elastic stress- 
intensity factor to "correct" for plastic yielding at the crack tip. The approach was to add 
the plastic-zone radius to the crack length and, thus, calculated a "plastic" stress-intensity 
factor at the "effective" crack length (c + r y ). The size of the plastic zone was estimated 
from equation (1) as 

r y = a i (Kfoys) (4) 

where ai = 11 ( 2 %) for plane stress and ai = l/(67t) for plane strain. The term aj is Irwin’s 
constraint parameter that accounts for three-dimensional stress state effects on yielding. 
Note that equation (4) gives approximately the "radius" of the plastic zone because of a 
redistribution of local crack-tip stresses due to yielding, which is not accounted for in the 
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elastic analysis. The actual plastic-zone size is roughly 2r y . Based on test experience, the 
fracture toughness, K c , calculated at the effective crack length remained nearly constant as 
a function of crack length and specimen width for several materials until the net-section 
stress exceeded 0.8 times the yield stress (a ys ) of the material [57]. 

Many researchers have used the Dugdale-Barenblatt (DB) model [55, 59] to 
develop some nonlinear crack-tip parameters (see Refs. 60 and 61). Drucker and Rice 
[62] presented some very interesting observations about the model. In a detailed study of 
the stress field in the elastic region of the model under small-scale yielding conditions, they 
reported that the model violates neither the Tresca nor von Mises yield criteria. They also 
found that for two-dimensional, plane-stress, perfect-plasticity theory, the DB model 
satisfies the plastic flow rules for a Tresca material. Thus, the model represents an exact 
two-dimensional plane-stress solution for a Tresca material even up to the plastic-collapse 
load. Therefore, the J-integral calculations from Rice [9] and AJ estimates may be 
reasonable and accurate under certain conditions. Of course, the application of the DB 
model to strain-hardening materials and to plane-strain conditions may raise serious 
questions because plane-strain yielding behavior is vastly different than that depicted by 
the model. 

Rice [9] evaluated the J integral from the DB model for a crack in an infinite body 
and found that 


J = a G 5 = 8 ct 0 2 c/(7tE) £n[sec(7cS/2a 0 )] (5) 

where a 0 is the flow stress, 6 is the crack-tip-opening displacement, c is the crack length, 
E is the elastic modulus, and S is the applied stress. An equivalent plastic stress-intensity 
factor Kj is given as 

Kj 2 = JE/(1 - n 2 ) (6) 

where r\ = 0 for plane stress, and r\ = v (Poisson's ratio) for plane strain. DB model 
solutions for plastic-zone size, p, and crack-tip opening displacement, 6, are available for a 
large number of crack configurations (see Ref. 41). Thus, J and Kj can be calculated for 
many crack configurations. However, for complex crack configurations, such as a 
through crack or surface crack at a hole, closed-form solutions are more difficult to 
obtain. A simple method was needed to estimate J for complex crack configurations. A-- 
common practice in elastic-plastic fracture mechanics has been to add a portion of the 
plastic zone p to the crack length, like Irwin's plastic-zone correction [56], to approximate 
the influence of crack-tip yielding on the crack-driving parameter. 

Newman [63] defined a plastic-zone corrected stress-intensity factor as 

K p = S (icd) 1/2 F(d/w,d/r, ...) (7) 

where d = c + yp and F is the boundary-correction factor evaluated at the effective crack 
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length. The term y was assumed to be constant and was evaluated for several two- and 
three-dimensional crack configurations by equating K p to Kj. From these evaluations, a 
value of 1/4 was found to give good agreement up to large values of applied stress to flow 
stress ratios. To put the value of one-quarter in perspective, Irwin's plastic-zone corrected 
stress-intensity factor [56] is given by y equal to about 0.4 and Barenblatt's cohesive 
modulus [59] is given by y = 1. A comparison of (elastic stress-intensity factor) and 
Kp, normalized by Kj and plotted against S/o G (applied stress to flow stress) for two 
symmetrical through cracks emanating from a circular hole is shown in Fig. 7. The dashed 
curves show Kg and the solid curves show K p for various crack-length-to-hole-radii (c/r). 
The elastic curves show significant deviations while the results from the K p equation (eqn. 
7) are within about 5% of Kj up to an applied stress level of about 80% of the flow stress. 


Through cracks at circular hole 



s/0 o 


FIG. 7— Ratio of elastic Kq and plastic Kp values to equivalent Kj values (after Newman, 
1992 [63]). 

This matches well with the 80%-limit established for Irwin’s plastic-zone corrected stress- 
intensity factor, as discussed by McClintock and Irwin [57]. To convert K p to AK p in 
equation 7, the applied stress and flow stress are replaced by AS and 2 g 0 , respectively, 
and p is replaced by the cyclic plastic zone co (see Ref. 64). Thus, Fig. 7 would be 
identical for cyclic behavior if Kp/Kj, is replaced by AKp/AKj and S/a 0 is replaced by 
AS/(2a 0 ), again, with y = 0.25. Thus, AK p is evaluated at a crack length plus one-quarter 
of the cyclic plastic zone. An application using this parameter to predict the fatigue life 
under high stresses will be presented later. 
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NUMERICAL ANALYSES AND MODELS OF CRACK GROWTH AND CLOSURE 


Since the early 1970’s, numerous finite-element and finite-difference analyses have been 
conducted to simulate fatigue-crack growth and closure. These analyses were conducted to 
obtain a basic understanding of the crack-growth and closure processes. Parallel to these 
numerical analyses, simple and more complex models of the fatigue-crack growth process were 
developed. Although the vast majority of these analyses and models were based on plasticity- 
induced crack-closure phenomenon, a few attempts have been made to model the roughness- 
and oxide-induced crack-closure behavior (see for example Refs. 36 and 65). This section will 
briefly review: (1) finite-element and finite-difference analyses, (2) yield-zone and empirical 
crack-closure models, and (3) the modified Dugdale or strip-yield models. In each category, an 
example of the results will be given. 

Finite-Element and Finite-Difference Analyses 

A chronological list of the finite-element and finite-difference analyses [66-88] is given in 
Table 1 . The vast majority of these analyses were conducted using two-dimensional analyses 
under either plane-stress or plane-strain conditions. Since the mid-1980’s, only a few three- 
dimensional finite-element analyses have been conducted. Newman and Armen [66-68] and 
Ohji et al. [69] were the first to conduct two-dimensional finite-element analyses of the crack- 
closure process. Their results under plane-stress conditions were in quantitative agreement 
with the experimental results of Elber [77] and showed that crack-opening stresses were a 
function of R ratio (S min/ S ma y) and stress level (S max /a 0 ). Nakagaki and Atluri [70] 
conducted crack-growth analyses under mixed-mode loading and found that cracks did not 
close under pure Mode II loading. In the mid- 1 980’ s, there was a widespread discussion on 
whether fatigue cracks would close under plane-strain conditions, i.e. where did the extra 
material to cause closure come from in the crack-tip region, since the material could not 
deform in the thickness direction, like that under plane-stress conditions. Blom and Holm [72] 
and Fleck and Newman [76, 79] studied crack-growth and closure under plane-strain 
conditions and found that cracks did close but the crack-opening levels were much lower than 
those under plane-stress conditions. Sehitoglu and his coworkers [74, 55] found later that the 
residual plastic deformations that cause closure came from the flanks of the crack (i.e., the 
material was plastically stretched in the direction parallel to the crack surfaces). Nicholas et al. 
[77] studied the closure behavior of short cracks and found that at negative R ratios the crack*" 
opening levels were negative, i.e. the short cracks were open at a negative load. 

In 1992, Llorca [84] was the first to analyze the roughness-induced closure mechanism 
using the finite-difference method. He found that the key controlling factor in roughness- 
induced closure was the tilt angle (0) between the crack branches (as the crack zig-zags ±0 
degrees). Crack-opening loads as high as 70% of the maximum load were calculated and these 
results agree with the very high opening loads measured on the 2124-T35 1 aluminum alloy. 
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TABLE 1-- Finite-element analyses of fatigue crack growth and closure . 


• Two-Dimensional Cracks 


• Three-Dimensional Cracks 


Newman 1974-76 Chermahini 1986 

Newman and Armen 1975 Chermahini et al. 1988-93 

Ohji, Ogura and Ohkubo 1975-77 Dawicke et al. 1992 

Nakagaki and Atluri 1980 Newman 1993 

Anquez 1983 

Blom and Holm 1985 

Kobayashi and Nakamura 1985 

Lalor, Sehitoglu and McClung 1986-87 

Fleck 1986 

Bednarz 1986 

Nicholas et al. 1988 

Fleck and Newman 1988 

Llorca and Galvez (a) 1989 

Anquez and Baudin 1988 

McClung et al. 1989-94 

Llorca (a) 1992 

Sehitoglu et al. 1993-96 


(a) Finite-difference method of analysis. 


McClung [80-82] performed extensive finite-element crack-closure calculations on small 
cracks, crack at holes, and various fatigue-crack growth specimens. Whereas Newman [65] 
found that Smax/ao could correlate the crack-opening stresses for different flow stresses (a 0 ) 
for the middle-crack tension specimen, McClung found that K-analogy, using Kmav/Kn could 
correlate the crack-opening stresses for different crack configurations for small-scale yielding 
conditions. The term Kq = g 0 V(7ia) where o 0 is the flow stress. (K-analogy assumes that the 
stress-intensity factor controls the development of closure and crack-opening stresses, and that 
by matching the K solution among different cracked specimens, an estimate can be made for 
the crack-opening stresses.) Some typical results are shown in Fig. 8. The calculated crack- 
opening stress, S op /S max , is plotted against K max /Ko for three crack configurations: middle- 
crack tension M(T), single-edge-crack tension SE(T) and bend SE(B). The symbols show the 
finite-element calculations for three crack-length-to-width ratios (c/w). At high values of 
Kn^/Ko, the crack-opening values became a function of crack configuration. A similar 
approach using "plastic-zone analogy" may be able to correlate the crack-opening stresses 
under large-scale yielding. The dashed curve shows the crack-opening stress equation, from 
the strip-yield model, developed by Newman [89] and recast in term of K-analogy. The 
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dashed curve gives a lower bound to the finite-element results. The reason for the differences 
between the finite-element and strip-yield model results must await further study. 
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FIG. 8— Configuration effects on crack-opening stresses (after McClung, 1994 [ 82 ]). 


Very little research on three-dimensional finite-element analyses of crack closure has 
been conducted. In 1986, Chermahini et al. [ 86 - 88 ] was the first to investigate the three- 
dimensional nature of crack growth and closure. He found that the crack-opening stresses 
were higher near the free surface (plane stress) region than in the interior, as expected. Later, 
Dawicke et al. [ 90 ] obtained experimental crack-opening stresses, similar to Chermahini’ s 
calculations, along the crack front using Sunder’ s striation method [ 91 ], backface-strain gages, 
and finite-element calculations. 

In reviewing the many papers on finite-element analyses, a few analysts were extending 
the crack at "minimum" load, instead of at maximum load for various reasons. Real cracks do 
not extend at minimum load and crack-closure and crack-opening behavior calculated from 
these analyses should be viewed with caution. As is obvious from Table 1, further study is 
needed in the area of three-dimensional finite-element analyses of crack growth and closure to 
rationalize the three-dimensional nature of closure with respect to experimental measurements 
that are being made using crack-mouth and backface-strain gages. Because the measurements 
give a single value of crack-opening load, what is the relation between this measurement and 
the opening behavior along the crack front? Is the measurement giving the free surface value, 
i.e. the last region to open? The crack-opening value in the interior is probably the controlling 
value because it is dominant over a large region of the crack front [ 87 , 90 ]. Also, more 
analyses are needed on the other forms of closure, such as roughness-induced closure. From 
the author’s point of view, plasticity- and roughness-induced closure work together to close 
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the crack and the phenomena are difficult to separate. In Llorca’s analyses [84] in the near- 
threshold regime, the plastic-zone size was smaller than the mesh points in the finite-difference 
method. Is the method able to accurately account for the mixed-mode deformation under these 
conditions? Residual plastic deformations in the normal and shear directions are what causes 
the crack surfaces to prematurely contact during cyclic loading. 

Yield-Zone and Empirical Crack-Closure Models 

A list of some of the more popular yield-zone models [92-97] and empirical crack- 
closure models [98-102] is given in Table 2. The Wheeler [92] and Willenborg et al. [93] 
models were the first models proposed to explain crack-growth retardation after overloads. 
These models assume that retardation exists as long as the current crack-tip plastic zone is 
enclosed within the overload plastic zone. The physical basis for these models, however, is 
weak because they do not account for crack-growth acceleration due to underloads or 
immediately following an overload. Chang and Hudson [705] clearly demonstrated that 
retardation and acceleration are both necessary to have a reliable model. Later models by 
Gallagher [94], Chang [95] and Johnson [96] included functions to account for both 
retardation and acceleration. A new generation of models was introduced by Bell and 
Wolfinan [98], Schijve [99], de Koning [700], Baudin et al. [707] and Aliaga et al. [702] that 
were based on the crack-closure concept. The simplest model is the one proposed by Schijve, 
who assumed that the crack-opening stress remains constant during each flight in a flight-by- 
flight sequence. The other models developed empirical equations to account for retardation 
and acceleration, similar to the yield-zone models. 


TABLE 2— Empirical yield-zone or crack-closure models . 


• Yield-Zone Models 

Willenborg et al. 1971 
Wheeler 1972 
Gallagher (GWM) 1974 
Chang (EFFGRO) 1981 
Johnson (MPYZ) 1981 
Harter (MODGRO) 1988 


• Crack-Closure Models 

Bell and Wolfman 1976 
Schijve 1980 

de Koning (CORPUS) 1981 
Baudin et al. (ONERA) 1981 
Aliaga etal. (PREFAS) 1985 


Lazzeri et al. [104] conducted fatigue-crack growth tests on a middle-crack tension 
specimen under a flight-by-flight load history (ATR-spectrum) at a mean flight stress level 
(Sig) of 75 MPa. Tests results are shown in Fig. 9. These results show an initial high rate of 
growth followed by a slowing down of crack growth from 7 to 10 mm and then a steady rise in 
the overall growth rates until failure. This behavior is what Wanhill [105] calls "transient crack 
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growth" under spectrum loading. Lazzeri et al. then made comparisons of the predicted crack 
length against flights from four of the empirical models (CORPUS, PREFAS, ONERA, and 
MODGRO, see Table 2) and one strip-yield model (FASTRAN-II, to be discussed later). The 
predicted results are shown with symbols in Fig. 9. The MODGRO model was very 
conservative, while the other three empirical models gave essentially the same results but under 
predicted the flights to failure. The FASTRAN-II model predicted failure at about 15 % 
shorter than the test results, but this model came closer to modeling the "transient crack 
growth" behavior, as discussed by Wanhill. This behavior has been traced to the "constraint- 
loss" regime in thin-sheet materials by Newman [106\. 



FIG. 9—Comparison of predictions from various models on aircraft spectrum (after Lazzeri 
etal., 1995 [104]). 


Modified Dugdale or Strip- Yield Models 

A chronological list of the modified Dugdale or strip-yield models [ 107-125] is given in 
Table 3. Shortly after Elber [33] discovered crack closure, the research community began to 
develop analytical or numerical models to simulate fatigue-crack growth and closure. These 
models were designed to calculate the growth and closure behavior instead of assuming such 
behavior as in the empirical models. Seeger [. 107] and Newman [66] were the first to develop 
two types of models. Seeger modified the Dugdale model and Newman developed a ligament 
or strip-yield model. Later, a large group of similar models were also developed using the 
Dugdale model framework. Budiansky and Hutchinson [109] studied closure using an 
analytical model, while Dill and Saff [108], Fuhring and Seeger [111], and Newman [112] 
modified the Dugdale model. Some models used the analytical functions to model the plastic 
zone, while others divided the plastic zone into a number of elements. The model by Wang and 


17 



Blom [118] is a modification of Newman’s model [112] but their model was the first to include 
weight-functions to analyze other crack configurations. All of the other models in Table 3 are 
quite similar to those previously described. The models by Nakai et al. [113], Tanaka [11 6] 
and Sehitoglu et al. [55] began to address the effects of microstructure and crack-surface 
roughness on crack-closure behavior. 


TABLE 3— Modified Dugdale or strip-yield crack-closure models . 


Seeger 1973 

Newman 1974 

Dill and Saff 1976 

Budiansky and Hutchinson 1978 

Hardrath et al. 1978 

Fuhring and Seeger 1979 

Newman 1981, 1990 

Nakai etal. 1983 

Sehitoglu 1985 

Keyvanfar 1985 


Tanaka 1985 
Ibrahim 1986 

Wang and Blom 1987, 1991 
Chen and Nisitani 1988 
de Koning and Liefting 1988 
Keyvanfar and Nelson 1988 
Nakamura and Kobayashi 1988 
Daniewicz 1991 
ten Hoeve and de Koning 1995 
Sehitoglu et al. 1996 


Atypical modified Dugdale model is shown in Fig. 10. This model [110, 112] uses bar 
elements to model the plastic zone and the residual plastic deformations left as the crack grows. 
Three-dimensional constraint is accounted for by using the constraint factor, a. For plane- 
stress conditions, a is equal to unity and for plane-strain conditions, a is equal to 3. The 
constraint factor has been used to correlate constant-amplitude fatigue crack growth rate data, 
as will be discussed later. 

CONSTRAINT EFFECTS ON CRACK-GROWTH BEHAVIOR 

The importance of constraint effects in the failure analysis of cracked bodies has long 
been recognized by many investigators. Strain gradients that develop around a crack front 
cause the deformation in the local region to be constrained by the surrounding material. This 
constraint produces multiaxial stress states that influence fatigue-crack growth and fracture. 

The level of constraint depends upon the crack configuration and crack location relative to 
external boundaries, the material thickness, the type and magnitude of applied loading, and the 
material stress-strain properties. In the last few years, a concerted effort (see Refs. 126-128) 
has been undertaken to quantify the influence of constraint on fatigue-crack growth and 
fracture. To evaluate various constraint parameters, two- and three-dimensional stress analyses 
have been used to determine stress and deformation states for cracked bodies. The constraint 
parameters that are currently under investigation are (1) Normal stress, (2) Mean stress, (3) T- 
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FIG. 10— A typical Dugdale or strip-yield model for plasticity-induced closure (after Newman, 
1981 [. 112 ]). 


stress, and (4) Q-stress. In 1960, an elevation of the "normal" stress was used by Irwin [56] in 
developing equation (4) using only the K solution. This is similar to the constraint factor used 
in the modified strip-yield models (see Ref. 112). McClintock [129] and Rice and Tracey 
[ISO] considered the influence of the mean stress, a m = (ai +02 + C3)/3, on void growth to 
predict fracture. The mean stress parameter is currently being used in conjunction with three- 
dimensional (3D), elastic-plastic, finite-element analyses to characterize the local constraint at 
3D crack fronts, see for example Reference 131. 

In the early 1970’s, the fracture community realized that a single parameter, such as K or 
J, was not adequate in predicting the plastic-zone size and fracture over a wide range of crack 
lengths, specimen sizes, and loading conditions. At this point, "two-parameter" fracture 
mechanics was bom. The second parameter was "constraint. " The characterization of 
constraint, however, has been expressed in terms of the next term(s) in the series expansion of 
the elastic or elastic-plastic crack-tip stress fields. In 1975, Larsson and Carlsson [132] 
demonstrated that the second term, denoted as the T-stress (stress parallel to the crack 
surfaces), had a significant effect on the shape and size of the plastic zone. The effects of the 
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elastic T-stress on J dominance for an elastic-plastic material under plane-strain conditions was 
studied by Betegon and Hancock [753] using finite-element analyses. Analytically, Li and 
Wang [134] developed a procedure to determine the second term in the asymptotic expansion 
of the crack-tip stress field for a nonlinear material under plane-strain conditions. Similarly, 
O’Dowd and Shih [735, 136] have developed the J-Q field equations to characterize the 
difference between the HRR stress field and the actual stresses. The Q-stress collectively 
represents all of the higher order terms for nonlinear material behavior. The J-Q field equations 
have been developed for plane-strain conditions. An asymptotic analysis that includes more 
terms for the stress and deformation fields at a crack embedded in a nonlinear material under 
Mode I and II loading for either plane-stress or plane-strain conditions has been developed by 
Yang et al. [737, 138]. Chao et al. [73<5] demonstrated that the first "three" terms in these 
series expansion (J, A 2 , and A 3 ) can characterize the stress Oy for a large region around the 
crack tip for Mode I plane-strain conditions. The third term was subsequently shown to be 
directly related to the first and second terms, thus two amplitudes, J and A 2 , were sufficient to 
describe the local stress field. 

In 1973, the Two-Parameter Fracture Criterion (TPFC) of Newman [139, 140] was 
developed which used the additional term in the local stress equations for a sharp notch or a 
crack. The TPFC equation, Kp = Kj e / (1 - m Sn/aJ, was derived using two approaches. (Kp 
and m were the two fracture parameters; Kj e is the elastic stress-intensity at failure; S n is the 
net-section stress and c u is the ultimate tensile strength.) In the first approach, the stress- 
concentration factor for an ellipsoidal cavity, Kp = 1 + 2 (a/p) /O from Sadowsky and 
Sternberg [141], was used with Neuber’s equation [ 6 ], K e = Kf , to derive a relation 
between local elastic-plastic stresses and strains and remote loading. This is similar to the way 
Kuhn and Figge [142] used the Hardrath-Ohman equation [7] many years earlier. Assuming 
that fracture occurred when the notch-root stress and strain was equal to the fracture stress and 
strain, af and ef, respectively, and that a crack had a critical notch-root radius, p*, the TPFC 
equation was derived. The second parameter, m, came from the "unity" term in the stress- 
concentration equation. The second approach [140], used the elastic stress field equation for a 
crack (eqn.(l)) and Neuber’s equation to relate the elastic stresses to the elastic-plastic stresses 
and strains at a crack tip. In this approach, it was again assumed that fracture occurred at a 
critical distance, r*, in front of the crack tip when the local stress and strain was equal to the 
fracture stress and strain, af and sf. The second parameter, m, came from the next higher- 
order term in the stress-field expansion. The TPFC has been successfully applied to a large 
amount of fracture data on two- and three-dimensional crack configuration and materials. 

As pointed out by Merkle, in the Fifth Swedlow Lecture [143], "... estimation of 
constraint effects is best accomplished with three-dimensional analyses. " With this in mind, 
Newman et al. [144, 145] conducted 3D elastic-plastic, finite-element analyses on a cracked 
plate with a wide range in crack lengths, thicknesses, and widths for an elastic-perfectly-plastic 
material under tension and bending loads. Because the previously discussed crack-closure 
models require information about constraint (elevation of the normal stress around the crack 
tip), an average normal stress in the plastically-deformed material normalized by the flow stress 
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was evaluated from the 3D analyses. This normalized average stress was denoted as a global 
constraint factor, <Xg. Some typical results of ctg plotted against a normalized K are shown in 
Fig. 1 1 for a thin-sheet material. The symbols show the results from the analyses for various 
specimen sizes. The upper dashed lines show the results under plane-strain conditions. The 
global constraint factor was nearly a unique function of the applied K level. Some slight 
differences were observed near the plane-stress conditions (high K levels). These results show 
that the global constraint factor rapidly drops as the K level increases (plastic-zone size 
increases) and approaches a value near the plane-stress limit. The solid line is a simple fit to the 
results and shows that the constraint-loss regime may be defined by a unique set of K values 
under monotonic loading. On the basis of some results from cyclic loading and conjecture, the 
constraint-loss regime may be defined by a unique set of AKgfj values under cyclic loading. 

This point will be discussed later. 



FIG. 1 1— Constraint variations from three-dimensional finite-element analyses (after Newman 
et al., 1994 [ 145 ]). 

CRACK GROWTH RATE RELATIONS 

The number of fatigue-crack growth rate relations in the literature is enormous. But the 
first such relation was attributed to Head in 1953 [ 146 ]. After the Comet accidents [2], which 
were caused by fatigue cracks growing from windows in the fuselage, the search for a reliable 
crack-tip parameter and growth rate relation was underway. Table 4 gives a very small list of 
some crack-growth rate relations that have been proposed since the early 1960’s. This list is a 
summary of the major relations that are currently being used today in many damage-tolerance 
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life-prediction codes. In 1961, Paris et al. [4] made a major step in applying the stress-intensity 
factor range to fatigue-crack growth. Donaldson and Anderson [ 147] demonstrated how this 
new concept could be applied to aircraft components. Very quickly it was found that AK alone 
would not correlate fatigue-crack growth rate data for different stress ratios, R, and other 
equations were proposed. Of these, the Forman et al. [148] and Walker [149] equations are 
commonly used in many life-prediction codes. The next major step in understanding fatigue- 
crack growth came when Elber [11, 33, 150] discovered crack closure and proposed that the 
AKgff parameter should control crack growth. Prior to Elber’ s discovery, Tomkins [151] 
was using the BCS (Bilby, Cottrell and Swinden [ 152 ]) model to develop a local crack-tip 
displacement parameter for crack growth. After Rice [9] developed the J-integral, 
Dowling and Begley [53, 54], and others, began to explore the use of the AJ e ff parameter 
for fatigue-crack growth. Similarly, Ogura et al. [753] proposed to use the local cyclic 
hysteresis energy (W e ff). The relationship between AK, or any other parameter discussed 
here, plotted against crack-growth rate does not always fit the simple power laws that 
have been proposed. Miller and Gallagher [154] found that more accurate life predictions 
could be made if a table-lookup procedure was used. (The reason that the table-lookup 
procedure is more accurate will become apparent later.) A number of life prediction 
codes, such as NASA FLAGRO [755] and FASTRAN-II [156], have adopted this 
procedure. 


TABLE 4— Evolution of some typical crack-growth rate relations . 

• Paris, Gomez and Anderson (1961): dc/dN = C AK n 

• Paris and Erdogan (1963): dc/dN = C AK n (K max ) m 

• Forman, Kearney and Engle (1967): dc/dN = f(AK, R, K c ) 

• Tomkins (1968): dc/dN = f(ACTOD) 

• Elber (1970): dc/dN = C (AK e ff)” 

• Walker (1970): dc/dN = f(AK, R) 

• Dowling and Begley (1976): dc/dN = C (AJ e ff) n 

• Ogura etal. (1985): dc/dN = f(W eff ) 

• Miller and Gallagher (1981): Table-lookup procedure dc/dN = f(AK, R) 


LARGE CRACK GROWTH BEHAVIOR 

This next section will review some observations and present results on the effects of 
microstructure, environment, and loading on fatigue-crack growth rate behavior. 
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Microstructural Effects 


As previously mentioned, fatigue-crack growth rate relations do not necessarily fit a 
simple power law because of sharp transitions in the AK-rate curves. In 1982, Yoder et al. 
[757], using data generated by Bucci et al. (1980), began to explain these transitions in terms of 
microstructural barriers to slip-band transmission, as shown in Fig. 12. The transition to 
threshold, Ti, appeared to be controlled by dispersoid spacing (cyclic plastic zone was about 
the size of the mean free path between dispersoid particles). Similarly, for T 2 and T 3 , the cyclic 
plastic zone size appeared to correlate with the subgrain and grain size, respectively. Note that 
these tests were conducted at an R ratio of 0.33. Phillips (see Ref. 145 ) tested 7075-T6 at R 
ratios of-1, 0 and 0.5 and found that these transitions were at different AK levels for each R 
ratio, but that each transition occurred at nearly the same crack-growth rate. Wanhill (Annex 
A in Ref. 17 ) found similar transitions in 2024-T3 aluminum alloy and concluded that the 
transitions were controlled by the "effective cyclic plastic zones" based on AK^g* instead of 
the cyclic plastic zones computed from AK. Thus, the transitions in the two alloys appear to be 
controlled by AKeff in laboratory air because the transitions occur at about the same rate. 


dc/dN 
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FIG. 12— Microstructural control of fatigue-crack growth (after Yoder et al., 1982 [757]). 
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Environmental Effects 


Piascik and Gangloff [755] found that these transitions were affected by the environment 
in crack-growth rate tests on a 2090 aluminum-lithium alloy. Fig. 13 shows tests results in (1) 
moist air or water vapor, (2) NaCl solution, and (3) Oxygen, helium or vacuum. Test results in 
each category fell along a particular AK-rate relation. The results in moist air or water vapor 
show a similar characteristic as exhibited by the 7075 alloy in laboratory air (Fig. 12), that is, 
the sharp transitions at T i and T 2 . But tests under the salt solution, eliminated the T 2 
transition and moved the Tj transition to a different values of AK. Under the inert 
environments, the transitions did not develop. Piascik and Gangloff attributed these behaviors 
at low crack-growth rates to fracture mode changes from cracking on the {100} plane in the 
salt solution to slip-band cracking in the inert environments. 
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FIG. 13— Environmental fatigue-crack growth in aluminum-lithium alloy (after Piascik and 
Gangloff, 1993 [755]). 


Results from Petit and Henaff [759, 160 ], as shown in Fig. 14, demonstrate the intrinsic 
behavior of crack growth under the various fracture modes. Stage I, Stage II and Stage I-Like 
for a wide variety of materials and R ratios tested in high vacuum. When the Stage II crack- 
growth rate data is plotted against AKgjf normalized by the elastic modulus (E), all materials 
(aluminum alloys, aluminum-lithium alloys, steels and TA6V) fall along a unique relation. 
Similar results are shown for the Stage I and Stage I-Like behaviors. These results help explain 
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why the fracture mode changes can produce transitions in the AK-rate relations. These results 
also demonstrate how new, metallic materials can be developed to have improved damage- 
tolerance properties. With many engineering metals falling together on a AK^/E plot, one 
way to improve the material is to produce a material that would have a large amount of closure 
either due to plasticity, roughness, or some other closure mechanism. 
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FIG. 14— Intrinsic fatigue crack growth for various material at high vacuum (after Petit and 
Henaff, 1991-93 [159,160]). 


Loading Effects 

Large-Crack Threshold — Because many of the comparisons between the growth of 
small and large cracks have been made in the near-threshold regime for large cracks, it is "~ 
important to know whether the large-crack threshold is a material property or is caused by 
the load-reduction procedure. Several investigators have experimentally or numerically 
shown (see Ref. 161) that the stress-intensity factor threshold under load-reduction 
schemes can be partly explained by the crack-closure behavior. Some typical results on an 
aluminum alloy are shown in Fig. 15. Minakawa and McEvily [162] conducted a 
threshold test on a compact specimen and measured the crack-opening loads as the AK 
level approached AK^. The crack-opening loads were determined from a displacement 
gage at the crack mouth. For high AK levels, the P(/Pmax values ranged from 0. 15 to 
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0.35. The horizontal line is the calculated Po/P ma x ratio from Newman’s crack-closure 
model [112] under constant-amplitude loading with plane-strain conditions (a = 3). The 
calculated ratios agreed fairly well with the experimental values. As AK approached AK&, 
the Pc/Pmax ratio rapidly rose and the ratio was nearly unity at threshold. Thus, the rise in 
crack-opening load explains why the threshold developed. But what caused the rise in 
crack-opening loads? A number of suggestions have been advanced to explain this 
behavior. Among these are the mismatch of crack-surface features observed by Walker 
and Beevers [34] in a titanium alloy; the corrosion product formation on the crack 
surfaces, as observed by Paris et al. [37]; the variation in the mode of crack growth with 
stress-intensity factor level as reported by Minakawa and McEvily [162]\ and the plastic 
wake caused by the load-reduction procedure [161]. The usefulness of the large-crack 
threshold data for small-crack growth is in question, if the threshold development is 
caused by the activation of different fracture modes, such as roughness-induced closure, 
and these mechanisms are not activated for small cracks. If the large-crack threshold is 
affected by the load-reduction procedure, then the overall usefulness of AK^ for large- 
crack-growth behavior under variable-amplitude loading is also questionable. 



FIG. 15— Experimental crack-growth rates and opening levels near large-crack threshold 
(modified after Minakawa and McEvily, 1981 [7(57]). 


Transition from Tensile-to-Shear Mode Crack Growth — The crack-growth regime 
where a crack grows from flat (tensile fracture) to slant (shear fracture), as shown in Fig. 
16, is important to defining the constraint-loss regime from plane strain to plane stress. As 
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observed by Schijve [26], the end of transition from flat-to-slant crack growth appears to 
occur at the same fatigue crack-growth rate, independent of the stress ratio. Newman et 
al. [163] used this observation to control the constraint-loss regime in the analytical crack- 
closure model. Because the crack-closure concept is able to collapse crack-growth rate 
data onto nearly a single AK e ff-rate relation, Schijve [164] proposed that the effective 
stress-intensity factor should control the transition from flat-to-slant crack growth. To 
develop a simple estimate for the transitional region, Newman [ 106] proposed that the 
transition to complete slant crack growth occurs when the effective cyclic plastic-zone size 
calculated from AK^g* is a certain percentage of the sheet thickness. This relation is 

p - (AKeff)x / (<* o (8) 

where a 0 is the flow stress and B is the sheet thickness. Using transitional data from the 
literature, the transitional coefficient (p) is plotted against sheet thickness in Fig. 17. 
Although considerable scatter is evident in the data, the general trend is for p to be about 
0.5 for 1 to 6 mm-thick material. While the AKeff at the end of transition is a function of 
specimen thickness, Wilhem [165] suggested that the beginning of the shear-lip 
development for aluminum and titanium alloys may be independent of specimen thickness. 
This is reasonable, considering that the material at the free surface is in a state of plane 
stress, regardless of thickness. 


Fatigue surface 



FIG. 1 6-Flat-to-slant fatigue-crack growth in metallic materials. 
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FIG. 17— Controlling parameter for flat-to-slant fatigue-crack growth (after Newman, 
1992 [106]). 


Constant-Amplitude Loading — In 1969, Hudson [ 166] produced fatigue-crack- 
growth-rate data for 2024-T3 and 7075-T6 aluminum alloy sheet over a wide range of 
stress ratios (R = -1 to 0.8) and stress-intensity factor ranges. Later, Phillips [145, 167] 
generated crack-growth data in the near-threshold regime for the same alloys; and 
Dubensky [168] conducted tests at extremely high remote stress levels (0.6 to 1.0 times 
the yield stress of the material). These tests produced crack-growth-rate data over 8- 
orders of magnitude in rates! These types of tests and data are needed to obtain the 
baseline crack-growth-rate relations that are needed to predict crack growth under 
variable-amplitude and aircraft spectrum loading, as will be discussed later. 

Typical fatigue-crack growth rate data on 7075-T6 aluminum alloy sheet for various 
R ratios [20] and analyzed with Newman’s closure model equations [89] are shown in Fig. 
18. On the basis of AKgff , the data collapsed into a narrow band with several changes in 
slope (transitions) occurring at about the same growth rate. For these calculations, a 
constraint factor (a) of 1.8 was selected for rates less than 7E-4 mm/cycle and a equal to 
1.2 for rates greater than 7E-3 mm/cycle. The vertical dash line shows the calculation of 
(AKgf^T from equation (8) with p = 0.5. The T4 location (defined herein) shows a sharp 
transition in the constraint-loss regime. The solid line is the baseline relation. In the low 
crack-growth rate regime, the large-crack threshold data has been neglected. The baseline 
relation near the large-crack threshold is an estimate based on small-crack data [20]. 
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FIG. 1 8— Effective stress-intensity factor against crack-growth rate for an aluminum alloy 
(after Newman et al., 1994 [20]). 

Spectrum Loading — Wanhill [169, 170] conducted spectrum crack-growth tests on 
middle-crack tension specimens made of 2024-T3 Alclad material (B = 3.1 mm). Tests 
were conducted under the TWIST (transport wing spectrum) clipped at Level III with a 
mean flight stress of S m f= 70 MPa. Fig. 19 shows a comparison of test results and 
calculated results from Newman’s closure model [112] with the constraint-loss regime (a 
= 2 to 1) estimated from equation (8). The model used AKgff-rate data like that shown in 
Fig. 18, but for the 2024-T3 alloy. To illustrate why the constraint-loss regime is 
necessary, example calculations were made for constant constraint conditions of either a - 
1 or 2 (dashed curves). The model with a low constraint condition (a = 1) predicted 
much longer lives than the tests, whereas the model with the high constraint predicted "~ 
much shorter lives than the tests. Thus, the correct constraint-loss regime is required to 
predict fatigue-crack growth under aircraft spectrum loading in thin-sheet materials. 

SMALL CRACK GROWTH BEHAVIOR 

The observation that small or short fatigue cracks can: (1) grow more rapid than 
those predicted by linear-elastic fracture mechanics (LEFM) based on large-crack data, 
and (2) grow at AK levels well below the large-crack threshold, has attracted considerable 
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FIG. 19— Experimental and calculated crack length against flights for an aluminum alloy 
(after Newman, 1992 [706]). 


attention in the last two decades [72-27]. Some consensus is emerging on crack 
dimensions, mechanisms, and possible methods to correlate and to predict small-crack 
behavior. A useful classification of small cracks has been made by Ritchie and Lankford 
[777] and these are summarized in Table 5. Naturally-occurring (three-dimensional) small 
cracks, often approaching microstructural dimensions, are largely affected by crack shape 
(surface or comer cracks), enhanced crack-tip plastic strains due to micro-plasticity, local 
arrest at grain boundaries, and the lack of crack closure in the early stages of growth. 
Whereas, two-dimensional short cracks, about 100 jum or greater, are through-thickness 
cracks which have been created artificially be removing the wake of material from large 
through cracks. Their behavior appears to be controlled by the plastic-wake history left by 
the large-crack growth process and the crack-growth rates are averaged over many grains 
through the thickness. 

Over the last two decades, in the treatment of microstructurally-, mechanically-, and 
physically-small cracks, two basic approaches have emerged to explain the rapid growth 
and deceleration of small cracks when compared to large-crack growth behavior. The first 
is characterized by "grain-boundary" blocking and consideration of microstructural 
effects on small-crack growth rates (see for example Refs. 27 and 772). The second is a 
"continuum mechanics" approach accounting for the effects of material nonlinearity on 
the crack-tip driving force and crack-closure transients (see for example Refs. 161 and 
173 ). 
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TABLE 5 — Classes of small-fatigue cracks, dimensions, responsible mechanisms and 
potential crack-tip parameters (modified after Ritchie and Lankford. 1986 |777]Y 


Types of Small Cracks 

Dimension 

Mechanisms 

Parameters 

Microstructurally-small 

a<d g (a) 

2c < 5 to 10 dg 

Crack-tip shielding 
enhanced Asp 
Crack shape 

Probabilistic 

approach 

Mechanically-small 

a<r y < b > 

Excessive (active) 
plasticity 

AJ, ACTOD 
(AKp) 

Physically-small 

a < 1 mm 

Crack-tip shielding 
(crack closure) 

AKeff 

Chemically-small 

up to = 10 mm 

(c) 

Local crack tip 
environment 



(a) d g is critical microstructural dimension, such as grain size; a is surface-crack depth 
and 2c is surface-crack length 

(b) ry is plastic-zone size or plastic field of notch 

(c) critical size is a function of frequency and reaction kinetics 


The microstructural barrier model, developed by Miller and co-workers [75, 772], 
was conceived to separate regimes of "microstructurally-small'' cracks and "physically- 
small" cracks. The regime of microstructurally-small cracks (MSC) occurs when crack 
lengths are less than a dominant microstructural barrier, such as the grain size. Various 
researchers consider this regime to be synonymous with growth of a crack across a single 
grain or several grain diameters. For example, a crack may initiate at an inclusion particle 
on a grain boundary, propagate, slow down, and stop at the next grain boundary. With 
further cycling, or if the stress level is increased, this barrier can be overcome and the 
crack will propagate to the next barrier. Several different microstructural barriers to crack 
growth may exist in a single material because of material anisotropy and texture. The ,~ 
physically-small crack (PSC) regime is defined for crack lengths greater than the spacing 
of these dominant barriers. Miller [772] suggests that the complexities near micro- 
structural barriers in the MSC and PSC regimes hinder theoretical analyses of small-crack 
growth behavior based on LEFM parameters and he emphasizes the development of 
empirical equations, based on extensive test data, to determine constants in these relations. 
However, progress has been made in the analyses of cracks growing from inclusions (see 
for example Ref. 174 ) and interacting with grain boundaries [775, 116 ]. These analyses 
may be useful in developing the LEFM relations for cracks in complex microstructures. 
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As shown in Fig. 20, small-crack initiation and growth is a three-dimensional 
process with cracks in the depth, a, and length, c, directions interacting with the grain 
boundaries at different times in their cyclic history. Whereas, an observed crack in the 
length direction may have decelerated at or near a grain boundary, the crack depth may 
still be growing. As the crack grows in the depth direction, the rise in the crack-driving 
force at the c-location contributes to the crack penetrating that barrier. As the cracks 
become longer, the influence of grain boundaries become less as the crack front begins to 
average behavior over more grains. Small-crack growth deceleration may or may not 
occur depending upon the orientation of the adjacent grains [27]. A probabilistic analysis 
would be required to assess the influence of the variability of the grain structure on crack- 
growth rate properties. From an engineering standpoint, however, a weak-link or worst 
case scenario of grain orientation may provide a conservative estimate for the growth of 
small cracks through a complex microstructure. This is the basis for the continuum 
mechanics approach. 

It has been argued that the calculation of AK for a small crack growing from an 
inclusion could be in error (Schijve [775]). For example, if crack initiation occurs at a 
subsurface inclusion with subsequent breakthrough to the surface, a considerable elevation 


Grain 

Crack front 


* c ’ 




FIG. 20— Surface crack growth and an influence of grain boundaries. 
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in AK is possible over that calculated from surface observations. Although the use of AK 
to characterize the growth of small cracks has proved to be convenient, its universal 
application has been viewed with some skepticism. Despite the abpve qualifications, 
research work on the growth of naturally-initiated small cracks, notably by Lankford [21, 
17 6] and AGARD studies [17, 18], has demonstrated the usefulness of the AK concept. 

One of the leading continuum mechanics approaches to small-crack growth is that of 
Newman et al. [161, 165]. The crack-closure transient (or more correctly the lack of 
closure in the early stages of growth) has long been suspected as a leading reason for the 
small-crack effect. The Newman crack-closure model [112] has demonstrated the 
capability to model small-crack growth behavior in a wide variety of materials and loading 
conditions [161]. Difficulties still exist for large-scale plastic deformations at holes or 
notches but these are problems that can be treated with advanced continuum mechanics 
concepts. In the remaining sections, the application of the Newman model to predict or 
calculate fatigue life based solely on crack propagation will be reviewed. 

Earlier work by Pearson [12] on fatigue-crack initiation and growth of small cracks 
from inclusion particles in two aluminum alloys (BS L65 and DTD 5050) set the stage for 
the development of small-crack theory. His results are shown in Fig. 21(a), as the dashed 
curve, along with additional small-crack data from Lankford [21] on 7075-T6 aluminum 
alloy. Pearson concluded that cracks of about the size of the average grain size, grew 
several times faster than large cracks at nominally identical AK values. The open symbols 
and dash-dot curve show the large-crack data and the development of the large-crack 
threshold at about 3 to 4 MPa Vm. The solid symbols show the typical small-crack 
behavior with growth at AK levels as low as 1 .5 MPa Vm. The other solid curves show the 
behavior measured on other small cracks. Some general observations, by Lankford [21], 
were that the minimum in da/dN occurred when the crack depth, a, was about the minimum 
dimension of the pancake grain (subsurface grain boundary, as shown in Fig. 20) and that 
the magnitude of the lower rates was controlled by the degree of micro-plasticity in the 
next grain penetrated by the crack. If the next grain is oriented like the first, then no 
deceleration will occur, as indicated by the uppermost small-crack curves in Fig. 21(a). 

At this stage, it would be of interest to compare the test results from Pearson and 
Lankford with the small-crack growth predictions made from a continuum-mechanics 
model based on crack closure [156]. The AK c ff-rate relation used in the closure model for 
the 7075-T6 alloy was generated in Ref. 20 and the relation is shown in Fig. 21(b) as the 
dotted lines. These results were generated from large-crack data for rates greater than — 
about 2E-6 mm/cycle. The results are quite different from those shown for the Pearson- 
Lankford large-crack data in that the Ti transition (at about IE-5 mm/cycle) did not 
appear in their data. The lower section of the AKefj-rate relation (below 2E-6 mm/cycle) 
was estimated on the basis of small-crack data, also generated in Ref. 20. Because small 
cracks are assumed to be fully open on the first cycle, the AKgff-rate relation is the 
starting point for small-crack analysis. The results of an analysis of the test specimen used 
by Lankford is shown by the heavy solid curve. The initial defect was selected as a 10 pm 
radius semi-circular crack, so that the 2a dimension (on the surface) would be 20 pm. 
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FIG. 21(a)— Measured and calculated small-crack growth in an aluminum alloy 
(after Lankford, 1982 [27]). 

As the small crack grows, the closure level increases much faster than the AK level and a 
rapid decrease in rates is calculated. This rapid drop is a combination of the closure 
transient and the sharp change in slope at the Tj location (at about IE-6 mm/cycle). At 
about 30 pm, the crack-opening stresses have nearly stabilized and the effects of plasticity 
on the crack-driving force is quite small considering that the applied stress was 0.75 times 
the flow stress (see Fig. 6 in Ref. 63). The predicted small-crack results are in excellent 
agreement with Pearson’s data and agree with Lankford’s data which do not exhibit a 
grain-boundary influence. Interestingly, the small-crack analysis shows a single dip in the 
small-crack curve, similar to the "single" dip observed in some of Lankford’s small-crack 
data. Would the grain-boundary interaction always occur at the same crack length (40 
pm)? Why aren’t there other dips, or small indications of a dip, in the rate curve at 80, 
120 or 160 pm? Further study is needed to help resolve these issues, but the fatigue life, 
or at least a lower bound fatigue life, can be calculated from continuum-mechanics 
concepts. The following sections will review the use of Small-Crack Theory to predict 
fatigue life for notch specimens under various load histories. 
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FIG. 21(b)— Measured and predicted small-crack growth in an aluminum alloy. 


PREDICTION OF FATIGUE LIFE USING SMALL-CRACK THEORY 

During the last decade, several international research programs have been conducted by 
the AGARD Structures and Materials Panel [77, 75, 177 , 178 ], the National Aeronautics and 
Space Administration (NASA) [ 20 ] and the Chinese Aeronautical Establishment (CAE) [ 179 ] 
on small-crack effects in a wide variety of materials and loading conditions. A summary of the 
materials and loading conditions used to generate small-crack data on notched specimens is 
listed in Table 6. Most of these studies dealt with naturally-initiated cracks at notches but some 
studies [ 177 , 178 ] used small electrical-discharged-machined notches to initiate cracks. 

All materials listed in Table 6 were subjected to a wide range of stress ratios (R) under 
constant-amplitude loading. The 2024-T3 aluminum alloy [ 17 , 75] was subjected to *- 

FAL STAFF, Inverted FALSTAFF, Gaussian, Felix-28, and TWIST loading; and the 7075-T6 
alloy [ 18 , 20 , 179 ] was subjected to the Gaussian and Mini-TWIST load spectra. The Lc9CS 
(a clad equivalent of 7075-T6) [ 20 , 179 ] was subjected to Mini-TWIST loading. Aluminum- 
lithium alloy 2090 [75] was subjected to the FALSTAFF, Gaussian, TWIST, and Felix-28 load 
sequences. The 4340 steel [75] was subjected to only the Felix-28 helicopter load sequence. 
The three titanium alloys (Ti-6A1-4V, IMI-685 and Ti-17) [777, 775] were subjected to only 
the Cold Turbistan spectrum. Details on these standardized spectra may be obtained from the 
appropriate references. 
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TABLE 6— Materials and loading conditions used in various small-crack programs . 


Materials 

Loading Conditions (a) 

2024-T3 

Constant-amplitude 

7075-T6 

FALSTAFF 

Lc9CS (clad) 

Inverted FALSTAFF 

2090-T8E41 

Gaussian 

4340 

TWIST 

Ti-6A1-4V 

Mini-TWIST 

IMI-685 

Felix-28 

Ti-17 

Cold Turbistan 
Commercial transport 


(a) List of materials and loading conditions are not associated. 


Newman and many of his coworkers [17, 20, 163, 180] used continuum-mechanics 
concepts with initial defect sizes, like those which initiated cracks at inclusion particles, voids 
or slip-bands, and the effective stress-intensity factor range (corrected for plasticity and 
closure) to predict the fatigue lives for many of the materials listed in Table 6. A summary of 
the initial defect sizes measured at numerous initiation sites is listed in Table 7. The baseline 
crack-growth rate data for these materials were obtained from the large-crack data, ignoring 
the large-crack threshold, and using small-crack growth rates at the extremely low rates. 
Small-crack thresholds were estimated from the endurance limit for the various materials. In 
the following, some typical examples of small-crack theory application will be presented. 

Aluminum Alloy 2024 

Landers and Hardrath [181] conducted fatigue tests on 2024-T3 aluminum alloy sheet 
material with specimens containing a central hole. The results are shown in Fig. 22 as symbols. 
Both elastic and elastic-plastic fatigue-crack growth analyses were performed. The solid and 
dashed curves show predictions using large-crack growth rate data (ignoring the large-crack 
threshold) and an initial crack size based on an average inclusion-particle size that initiated 
cracks [17]. The large-crack growth rate properties are given in Ref. 163 for the elastic stress- 
intensity factor analysis. The crack-growth rate properties using the elastic-plastic effective 
stress-intensity factor analysis were obtained from a re-analysis of the large-crack data [182]. 
Both predictions agreed near the fatigue limit but differed substantially as the applied stress 
approached the flow stress (a 0 = 425 MPa). In these predictions, a AK-effective threshold for 
small cracks was selected as 1 .05 MPa Vm (see Ref 17). Using this threshold value and the 
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TABLE 7— Average material defect sizes and equivalent area crack sizes 
(surface defect unless otherwise noted!. 


Defect 

Material Half-length, 

pm 


Equivalent Area Semi- 
Circular Crack Radius, 
Jim 

2024-T3 

3 

12 

6 

7075-T6 

3 

9 

5.2 

Lc9CS (clad) 

77 (a) 

77 

77 

2090-T8E41 

3.5 

11 

6 

4340 

8 

13 

10 

Ti-6A1-4V (sheet) 

(b) 

(b) 

0.5 (c) 

Ti-6A1-4V (forging) 

(b) 

(b) 

15 (c) 


(a) Comer crack at clad layer (clad thickness about 60-70 jim). 

(b) No metallurgical examination of initiation sites was made. 

(c) Value selected to fit fatigue data. 



Nf, cycles 


FIG. 22— Measured and predict fatigue lives for a hole in an aluminum alloy under constant- 
amplitude loading (after Newman, 1992 [ 182 ] ). 
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initial defect size, the analysis were used to match the fatigue limits for the two R ratio tests. 
Above a stress level of about 250 MPa, the results from the elastic and elastic-plastic analyses 
differed substantially. The results from the elastic-plastic analyses agreed well with the test 
data and substantiated the use of the cyclic-plastic-zone corrected effective-stress-intensity- 
factor range [ 63 ]. Similar comparisons between measured and predicted S-N behavior for 
notched and un-notched aluminum alloy specimens subjected to either constant- or variable- 
amplitude loading are presented in Reference 183 . 

Laz and Hillbeny [ 184 ] conducted a study to address the influence of inclusions on 
fatigue-crack formation in 2024-T3 using a probabilistic model in conjunction with the 
FASTRAN-H code [ 156 ]. An examination of the microstructure produced data on nearly 
3800 inclusions on the primary crack plane. A Monto Carlo simulation was conducted for 
1000 trials. In each trial, an inclusion particle area is randomly selected and converted to an 
equivalent semi-circular initial crack size. The initial crack size, which was assumed to be 
located at the center of a single-edge-notch tension (SENT) specimen, was grown using the 
FASTRAN-II model. The probabilistic model produced a predicted distribution of fatigue 
lives. These results are presented in Fig. 23 as the solid curve and are compared with 
experimental data from the AGARD study [ 17 ] and other tests conducted by Laz and 
Hillberry. The probabilistic model accurately predicts both the mean and variation of the 
experimental results. Most importantly, however, the model predicted the critical lowest life 
values, which correspond to the largest, most damaging inclusions. 



FIG. 23 —Measured and predicted cumulative distribution function for a notched aluminum 
alloy under constant-amplitude loading (after Laz and HUlbeny, 1995 [ 184 ]). 
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High-Strength 4340 Steel 


Swain et al. [755] conducted fatigue and small-crack tests on 4340 steel, single-edge- 
notch tensile specimens. Tests were also conducted on large cracks to obtain the baseline 
crack-growth rate data. These tests were conducted under both constant-amplitude and 
spectrum loading. Only the results from the spectrum load fatigue tests will be discussed here. 

The results of the fatigue life tests under the Felix-28 load sequence are shown in 
Fig. 24 as symbols. Inspection of the fracture surfaces showed that in each case a crack 
had initiated at an inclusion particle defect. Those specimens which had the shorter lives 
had cracks which initiated at spherical calcium-aluminate particle defects, whereas those 
with the longer lives had cracks which initiated at manganese sulfide stringer inclusion 
particle defects. Examination of the initiation sites for over 30 fatigue cracks produced 
information on the distribution of crack initiation site dimensions. The spherical particle 
defects range in size from 10 to 40 pm in diameter. The stringer particles were typically 
5 to 20 pm in the thickness direction and range up to 60 pm in length. The median values 
of the defect dimensions measured were aj = 8 pm and q = 13 pm. 



Nf, cycles 


FIG. 24—Measured and predict fatigue lives for notched high-strength steel under helicopter 
spectrum loading, Felix/28 (after Swain et al., 1990 [755]). 


Predictions of total fatigue life under the Felix-28 load spectrum were made using 
FASTRAN-II [756] by calculating the number of cycles necessary to grow a crack from 
the assumed initial defect size, located at the center of the notch root, to failure. The 
comparison between model predictions and experimental fatigue lives are shown in Fig. 
24. The predicted results agreed well for those specimens which contained spherical 
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defects as crack initiation sites. But, the predicted lives fell short for those specimens 
where cracks initiated from stringer inclusions or where no cracks of minimum size 
necessary for continued propagation were formed. Again, it is desirable to have good 
agreement for the larger spherical-inclusion particles because these particles will produce 
the minimum fatigue lives. For an engineering component, which may contain a large 
number of fastener holes or other areas of stress concentration, the likelihood of a critical 
sized inclusion particle being located at one of these sites is large. 

Titanium Alloy Ti-6A1-4V 

The titanium alloys listed in Table 7, the sheet and forging products, exhibited quite 
different behavior in terms of the equivalent initial flaw size (EIFS) required to fit the S-N data 
on notched specimens. For the sheet alloy, the EIFS was 0.5 pm [ 180 ], whereas the forging 
alloy required a defect 2 to 20 pm in depth to bound the scatter in tests on two different 
forgings. A comparison of the measured and calculated fatigue lives on the two titanium 
forgings are shown in Fig. 25. The fatigue test data generated on double-edge-notch tensile 
specimens made from the two forgings agreed well with each other. The FASTRAN-II life 
prediction code was used with large-crack growth rate data to calculate the fatigue lives from a 
selected initial semi-circular surface crack size, as shown by the solid curves (see Ref 145 ). An 
average defect size of about 10 pm would fit the mean of the experimental data quite well. 
(Note: An error was detected in the review process. The predicted S-N behavior. Fig. 20 in 
Ref. 145 , was calculated at R = 0 instead of 0. 1. Fig. 25 is a re-analysis of the S-N behavior.) 

Because no metallurgical examination of the fractured specimens was made to evaluate 
the initiation site in either References 178 or 186 , an attempt was made to relate the assumed 
initial crack size to some microstructural features (see Ref. 180 ). Because titanium has a 
relatively high solubility for most common elements and when multiple vacuum arc melting is 
accomplished with high purity materials, the occurrence of inclusion-type defects is rare. In a 
mill-annealed titanium alloy, Eylon and Pierce [ 187 ] studied the initiation of cracks and found 
that cracks preferred to nucleate in the width direction of alpha needles, or colonies of alpha 
needles, along a shear band on the basal plane. The size of the alpha needles was not reported 
but was stated to be considerably smaller than the size of alpha grains (grains were about 8 pm 
in size). Thus, the alpha needle size may be close to the initial crack size needed to predict 
most of the fatigue life based solely on crack propagation. For the forging alloy, Wanhill and 
Looije [755] found that the primary a-grains and platelet a-packets ranged from 8 to 18 pm. 
Whether these a-grains or packets contributed to the early initiation of micro-crack growth in 
the titanium forgings must await further study. 

DESIGN CONCEPT USING SMALL-CRACK THEORY 

The research that has occurred over the last two decades on "small- or micro-" crack 
growth has evolved into a new design concept that provides an alternative to the traditional 
safe-life (or S-N) approach. The diagram in Fig. 26 shows the various design concepts: Safe- 


40 




Nf, cycles 


FIG. 25— Measured and calculated fatigue lives for a double-edge-notch tensile specimen under 
constant-amplitude loading (modified after Newman et al., 1994 [ 145 \). 


life. Small-crack theory. Durability, Damage tolerance, and Fail-safe, and their relative location 
with respect to flaw size. The flow in the diagram, from safe-life to the fail-safe concept, 
depicts crack formation, micro-crack growth, macro-crack growth, and fracture. 

One of the original design concepts used in the aerospace industry was "safe-life", see 
Reference 22. Here a structure was assumed to be defect free and the life was established by 
conventional S-N or e-N approaches. A safe-life component is retired from service when its 
useful life has been expended. More recently, the total life has been calculated from a crack- 
initiation period (Nj) and a crack-propagation period (N p ), see Reference 189 . However, the 
crack size existing after the initiation period was debatable. In practice, however, it has been 
found that preexisting manufacturing defects (i.e., scratches, flaws, burrs, and cracks) or 
service induced damage (i.e., corrosion pits) are very often the source of structural cracking 
problems. The effect of these defects on the life of a component was dependent on the defects 
initial size, the rate of crack growth with service usage, and the critical crack size. 

From these considerations, the "damage-tolerance" design philosophy [ 190 ] was 
developed. The damage-tolerance approach assumes that a component has prior damage or a 
flaw size that is inspectable. Useful service life is then some specified portion of the crack 
propagation life required to grow the crack from the inspectable size to failure. 

The durability design concept [ 191 ] is economic-life extension to minimize in-service 
maintenance costs and to maximize performance. The durability analysis methodology, based 
on a probabilistic fracture mechanics approach, accounts for the initial fatigue quality (i.e., 
surface treatment, manufacturing defects), fatigue-crack growth in a population of structural 
details, load spectra, and structural design details. A statistical distribution of the "equivalent 
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FIG. 26—Evolution of design concepts in relation to "small-crack theory. " 


initial flaw sizes" (EIFS) is used to represent the initial fatigue quality of structural details. In 
practice, the EIFS distributions have ranged from 0. 1 to 0.5 mm. As the EEFS’s become 
smaller (below about 1 mm), the question of small-crack effects is raised. The current 
durability methodology, which does not account for small-crack effects, masks the small-crack 
behavior by deriving an EIFS that predicts the desired crack-propagation life. However, EIFS 
variability with load spectra may be resolved if small-crack theory is implemented into the 
design practice. 

The evolution of "safe-life" and "durability/damage-tolerance" design concepts have 
both moved in the direction of "small-crack behavior." As inspection techniques and 
manufacturing quality improves, smaller flaw sizes will be detected or produced in 
manufacturing. The design of uninspectable or hidden structures, or structures subjected to an 
extremely large number of cyclic loads (i.e., engines and helicopter components) must relie on 
approaches that deal with the growth of small cracks. 

"Small-crack theory" is currently used to assess the structural fatigue life at two levels 
of initial flaw sizes. First, it is used to assess the initial design quality of a structure based solely 
on "material" microstructural properties, such as cracks growing from inclusion particles, 
voids, grains, grain boundaries or cladding layers. Of course, this crack-growth life is the best 
the current material can provide under the desired loading. In the second level of calculations, 
as in the durability analysis, a manufacturing or service-induced flaw size is used with small- 
crack theory. But the effects of small-crack growth are now accounted for in the analysis (i.e., 
faster growth rates at a given AK and growth below the large-crack threshold). Small-crack 
theory can now be used to evaluate EIFS for various spectra and structural details; and to 
characterize S-N or e-N behavior. The impact of small-crack effects on design-life calculations 
have been discussed by Phillips and Newman [ 192 ]. 
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SUMMARY — Past, Present and Future 


A summary of the observations as well as past, present and potential future approaches 
concerning fatigue concepts, crack-propagation concepts, and small- (or micro-) crack growth 
behavior has been made. In the past, fatigue was characterized by stress-life (S-N) or strain-life 
(e-N) curves; and crack propagation was characterized by the AK-rate concept. It was 
observed that small cracks initiated early in life at high stress levels, but initiated late in life at 
low stress levels (near the endurance limit). Thus, initiation life was dominate near the 
endurance limit. At present, fatigue is characterized by a crack initiation stage (Nj) and a 
crack-propagation stage (N p ); and crack propagation is characterized by the effective cyclic 
stress-intensity factor (AKeg) or J-integral (AJ e ft). Small cracks grow faster than large cracks; 
and small cracks grow at stress-intensity factor (AK) levels well below the large-crack 
threshold (AKth)- 

Proposed future approaches are: (1) fatigue damage is characterized by crack size and 
fatigue lives are calculated by crack propagation from a micro-crack size, (2) crack growth is 
characterized by a local nonlinear crack-tip stress or deformation parameter, such as J, T*, 
ACTOD or cyclic hysteresis energy, (3) small-crack growth is viewed as the typical behavior 
and large-crack growth near and at threshold is the anomaly on a AK basis, and (4) fatigue life 
prediction methods need to be developed to predict micro-crack growth, as influenced by 
microstructure and environment, under complex load histories. The proposed approaches will 
potentially bridge the gap between "safe life" and "durability/damage tolerance" concepts. 
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mechanics" concepts; and has bridged the gap between safe-life and durability/damage-tolerance design concepts. 
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